Entropy-induced separation of star polymers in porous media 
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We present a quantitative picture of the separation of star polymers in a solution where part 
of the volume is influenced by a porous medium. To this end, we study the impact of long-range- 
correlated quenched disorder on the entropy and scaling properties of /-arm star polymers in a good 
solvent. We assume that the disorder is correlated on the polymer length scale with a power-law 
decay of the pair correlation function g(r) ~ r~ a . Applying the field-theoretical renormalization 
group approach we show in a double expansion in e = 4 — d and 8 — 4 — a that there is a range of 
correlation strengths S for which the disorder changes the scaling behavior of star polymers. In a 
second approach we calculate for fixed space dimension d — 3 and different values of the correlation 
parameter a the corresponding scaling exponents 7/ that govern entropic effects. We find that 
7/ — 1, the deviation of 7/ from its mean field value is amplified by the disorder once we increase 
8 beyond a threshold. The consequences for a solution of diluted chain and star polymers of equal 
molecular weight inside a porous medium are: star polymers exert a higher osmotic pressure than 
chain polymers and in general higher branched star polymers are expelled more strongly from the 
correlated porous medium. Surprisingly, polymer chains will prefer a stronger correlated medium to 
a less or uncorrelated medium of the same density while the opposite is the case for star polymers. 

PACS numbers: 64.60.Fr,61.41.+e,64.60.Ak,11.10.Gh 



I. INTRODUCTION 

The influence of structural disorder on the scaling 
properties of polymer macromolecules, dissolved in a 
good solvent is subject to ongoing intensive discussions 

E S S S S i II I S 113 IB El For polymers, 
structural disorder may be realized experimentally by a 
porous medium. Depending on the way the latter is pre- 
pared, it can mimic various behavior, ranging from un- 
correlated defects 0, 0, 0, to complicated fractal 
objects 0,EHH3- Consequently, theoretical and Monte 
Carlo (MC) studies have considered these different types 
of disorder. In particular, the scaling properties of poly- 
mer chains were analyzed for the situations of weak un- 
correlated 0, 3 ) of long-range-correlated 0, 0, as 
well as of fractal disorder at the percolation threshold 

I! ill SI EH 

However, as far as the authors know, 
the influence of correlated disordered media on the be- 
havior of branched polymers, e.g. polymer stars, have 
found less attention. Our work is intended to fill this 
gap. 

The study of star polymers is of great interest since 
it has a close relationship to the subject of micellar and 
other polymeric surfactant systems pll I22I l23| . More- 
over, it can be shown, that the scaling behavior of sim- 



ple star polymers also determines the behavior of general 
polymer networks of more complicated structure |24ll2flj . 
Recently, progress in the synthesis of high qu ality mono- 
disperse polymer networks |2fj, 1271 l28l l29l |3fJ has stimu- 
lated numerous theoretical studies of star poly mers, both 
by computer simulation [sil I32I Is!! Is^. Is^lj ^fjl] an d by the 
renormalization group techniqu e l24i l25l 1371 IsH Isil l4?iL 
EllIllSHISlSlillii, HIT Let us note, that 
polymer stars are hybrids between polymer like entities 
and colloidal particles, establishin g an im por tant link be- 
tween these different systems PII2II2II50II5II52I. 
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FIG. 1: Polymer star with / arms. 

It is well established that long flexible polymer chains 
in good solvents display universal and self-similar confor- 
mational properties on a coarse-grained scale and that 
these are perfectly described within a model of self- 
avoiding walks (SAWs) on a regular lattice |U 13 EH • 
For the average square end-to-end distance R and the 
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number of configurations Zm of a SAW of N steps one 
finds in the asymptotic limit N — > oo: 

(R 2 ) ~ A 2ly , Zjv ~ e^ATf" 1 (1) 

where and 7 are the universal exponents depending 
only on the space dimensionality d, and is a non- 
universal fugacity. The universal properties of this poly- 
mer model can be described quantitatively with high pre- 
cision by analyzing a correspon ding field the ory by renor- 
malization group method s IH4I l55l l56l I5H l58l l53 | . For 
d = 3 the exponents read |60|^ (0) = 0.5882 ±0.0011 and 
7^ = 1.1596 ± 0.0020. Here, and in the following we use 
the notation for the value of an exponent derived for 
the pure solution without disorder. 

The power laws of Eq. can be generalized to de- 
scribe a star polymer that consists of / linear polymer 
chains or SAWs, linked together at their end-points (see 
Fig. 1). For a single star with / arms of N steps 
(monomers) each, the number of possible configurations 
scales according to |2J, |2£j : 

Z NJ ~ e^AH/-! _ ( R /iyi,-fn 2 ( 2 ) 

in the asymptotic limit N — > 00 . The second part shows 
the power law in terms of the size R ~ N v of the isolated 
chain of N monomers on some microscopic step length 
i, omitting the fugacity factor. The exponents 7/, r\f 
are universal star exponents, depending on the number 
of arms /. The relations between these exponents read 


7/ = 1 + "(vf - fm) 

7i=72=7 = l-vrfr, Vi = 0. (3) 

Here, v and 7 are usual SAW exponents (JJJ. For / = 1, 
2 the case of a single polymer chain is restored. Recent 
numerical values for Tffor different / at d = 3 are given 
in Refs. 0, HI HI HiUiyplfor Monte Carlo (MC) 
simulations and in Refej25|l|2|3 HI EJ El El El for 
renormalization group calculations. 

In terms of the mutual interaction, polymer stars inter- 
polate between single polymer chains (low /) and poly- 
meric micelles (high /) [HI 03, Is^ . From the scaling 
properties of star polymers, one may also derive their 
short distance effective interaction. The mean force F(r) 
between two star polymers of / and /' arms is inversely 
proport ional to the distance r < R between their cores, 



with the amplitude given by the universal contact ex- 
ponent 6//'- The contact exponents are related to the 
family of exponents r\f for single star polymers by the 
following scaling relation (0 24j : 

®ff>=Vf+Vf'-Vf+f>- (5) 




FIG. 2: Separation phenomena of polymer stars in good 
solution, part of which is in a porous medium. 

Similar to the model of SAWs on a regular lattice which 
is used to describe the scaling properties of long flexible 
polymer chains in a good solvent, one may consider mod- 
els of SAWs on disordered lattices to study polymers in a 
disordered medium. In this model, a given fraction of the 
lattice sites is randomly chosen to be forbidden for the 
SAW (these forbidden sites will be called defects here- 
after). Harris [3| conjectured that the presence of weak 
quenched uncorrelated point-like defects should not alter 
the SAW critical exponents. This was later confirmed 
by renormalization group considerations Another 
picture appears however for strong disorder, when the 
fraction of allowed sites is at the percolation threshold. 
Numerous data from exact enumeration, analytical and 
MC simulation H H 0, H H El strongly suggest that 
the scaling of a SAW on a percolation cluster belongs to 
a new universality class and is governed by exponents, 
that differ from those of a SAW on a regular lattice. 

Our present study concerns the scaling properties of 
star polymers in porous media which are found to display 
correlations on a mesoscopic scale . In small angle X- 
ray and neutron scattering experiments these correlations 
often express themselves by a power law behavior of the 
structure factor S(q) ~ q~ df on scales £ _1 < q < l~~ x 
where I is a microscopic length scale and £ is the corre- 
lation length of the material and df is its fractal volume 
dimension 63]. We describe this medium by a model of 
long-range-correlated (extended) quenched defects. This 
model was proposed in Ref. 64] in the context of mag- 
netic phase transitions. It considers defects, character- 
ized by a pair correlation function g(r), that decays with 
a distance r according to a power law: 

g(r) ~ r~ a (6) 

at large r. For the structure factor this leads to a power 
law behavior with fractal dimension df — d — a where d is 
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the Euclidean space dimension. This type of disorder has 
a direct interpretation for integer values of a. Namely, 
the case a = d corresponds to point-like defects, while 
a = d — 1 (a = d — 2) describes straight lines (planes) 
of impurities of random orientation. Non-integer values 
of a are interpreted in terms of impurities organized in 
fractal structures [64| . 

The influence of the long-range-correlated defects © 
on magnetic phase transitions has been pointed out in 
theoretical work El El EE El El E3] and MC sim- 
ulations [7(3, El u3- F° r polymers, its impact on the 
scaling of single polymer chains was analyzed in our pre- 
vious work in two complementary renormalization group 
approaches: first by a double expansion in the parame- 
ters e = 4 — d and the correlation strength 5 = 4 — a 
using a linear approximation [lOj and secondly by evalu- 
ating two-loop expressions of the theory for fixed values 
of a and d [13, • In particular, this work showed that 
long-range-correlated disorder leads to a new universality 
class with values of the polymer scaling exponents that 
depend on the strength of the correlation expressed by 
the parameters a or S = 4 — a. From this we may expect 
that also the architecture dependent scaling behavior @ 
of polymer stars and networks is affected by this type of 
correlated disorder. 

The question we are interested in is: how does the 
presence of long-range-correlated disorder change the val- 
ues of the critical exponents 10), ©? Besides the star- 
star interaction, the exponents govern various phenom- 
ena that involve star polymers and polymer networks 
E! |H 0E Eg. A particular effect that may 
be observable experimentally for star polymer solutions 
in a porous medium is an architecture-dependent impact 
of the medium on the star polymer. It may lead to a 
separation of star polymers with different numbers / of 
arms. Let us consider star polymers in a good solvent, 
part of which is in a porous medium (see Fig. 2). We 
consider the pores to be large enough, so that the star 
polymers may pass in and out of the medium (however, 

possibly on long time scales only). Be J^P\N) the free 
energy of a star polymer with / arms of N steps each 
in the pure solvent and its free energy in a porous 
medium characterized by a correlation strength S. These 
can be estimated using (0: 

T ( f\N) = -\nZf\N) 

= -p(0)Nf-( 1 f-l)lnN, (7) 
T ( f S) (N) = -\nZf(N) 

= -p(p)Nf-( 1 f-l)lnN. (8) 

Here, we assume the fugacity factor e~^ p ) to depend on 
the concentration of impurities independent of their cor- 
relation, as it would be the case for SAWs on a lattice 
with corresponding defects ||. The product Nf repre- 
sents the total number of steps or effective monomers of 
the star polymer which is a dimensionless measure of its 



molecular mass. Using Eqs. |0 and JSJ) one may now 
compare the free energies of a number of situations. Let 
us name mainly two specific questions: (i) Given a star 
polymer with fixed mass Nf and functionality / in a 
good solution in a volume that is influenced by disor- 
der with a fixed defect density p. Does the free energy 
depend on the correlation, and in particular is the un- 
corrected disorder or rather the correlated disorder of 
the same density favored by the star polymer? (ii) Given 
a mixture of star polymers which are mono-disperse in 
mass Nf but polydisperse in functionality / in a good 
solution in which only a part of the volume is influenced 
by defects (see Fig. 0. Due to the fugacity contribution 
which is the same for all these star polymers, they are 
expected to favor the pure part of the solution. However 
the extent to which this is the case may depend on ar- 
chitecture. Is the star polymer mixture partly separated 
in this situation and where is the concentration of higher 
branched star polymers enhanced in this case? While our 
answer to the first question is mainly to be compared with 
MC simulations of star polymers on disordered lattices, 
the answer to the second one may also be relevant for 
experiments with polymers in solutions inside correlated 
structures like aerogels. 

The setup of the paper is as follows. In the next section 
we present the model and construct the Lagrangean of 
the corresponding field theory. In section III we describe 
the field-theoretical renormalization group (RG) methods 
that we apply. Section IV presents our results for the two 
RG approaches. We conclude with an interpretation of 
these results in section V. 



II. THE MODEL 

Let us consider a single star polymer with / arms im- 
mersed in a good solvent (Fig. 1). Working within the 
Edwards continuous chain model |73, , we represent 
each arm of the star by a path r a (s), parameterized by 
< s < S, a = 1, 2, . . . , /. In a corresponding discrete 
model of chains with TV steps of mean square microscopic 
length I the so-called Gaussian surface is S = N£ 2 . The 
central branching point of the star is fixed at i"i(0). The 
partition function of the system is then defined by the 
functional integral [25j: 

Zf(S) = J D[T U ...,T f ] 

f 

xexp[-H f ] Ip d (r Q (0)-n(0)). (9) 

a=1 

Here, Hf is the Hamiltonian, describing the system of / 
disconnected polymer chains: 
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+i E / ds / ds'6 d (r a (s)-r a ,(s')). 

a,a ; — 1 ^ ^ 

The first term in llOfl represents the chain connectivity 
whereas the second term describes the short range ex- 
cluded volume interaction. The product of <5-functions 
in 10 ensures the star-like configuration of the set of / 
chains requiring each of them to start at the point i"i(0). 
This model may be mapped to a field theory by a Laplace 
transformation from the Gaussian surface S to the con- 



jugated chemical potential variable (mass) £iq: 

Z f (fio)= J dSexp[-fi Q S}Z f (S). (11) 

One may then show that the Hamiltonian TL is related 
to an m-componcnt field theory with a Lagrangean C in 
the limit m — > and that the partition function Zf(p,o) 
results from a correlation function of this field theory as 
follows: 



Zf(th) = I d d Xl ---d d x f ( T ll _ lf ^(x )---^(x a )^(x 1 )---^(x f ))i^ , (12) 



d d x [ {»l\${x)\ 2 + \X7$(x)\ 2 ) + ^S^.^ix) . . . 4h(x) 



(13) 



Here and below, the summation over repeated indices 
is implied, <f> is an m-component vector field <f> — 
(tfi , . . . , 4> m ), and uq are bare mass and coupling with 

the tensor Si 1 __ j i 4 — ^ (^1*2^324 ^1*3^2*4 ^1^4^2*3)* 
Formally, the local composite operator appearing in Eq. 
(|12[) is the m = limit of an operator known in m- 
component field theory |75| : 

[<f>]l(x)=f il _ if ^(x)...^(x) (14) 

where T^,.,^, is a traceless symmetric SO(m) tensor: 

m 

E 7 '- "• ( 15 ) 

i=l 

We introduce disorder into the model l(T3)) . by redefin- 
ing (3lq — > fiQ + Sflo(x), where the local fluctuations 6fio(x) 
obey: 

((Sfi (x)}) = 0, 

{(Sfio(x)dfio(y))) = g{\x - y\) . 

Here, ((•••)) denotes the average over spatially homoge- 
neous and isotropic quenched disorder. The form of the 
pair correlation function g(r) is chosen to decay with dis- 
tance according to the power law ©■ 

In order to average the free energy over different con- 
figurations of the quenched disorder we apply the replica 
method to construct an effective Lagrangean: 

1 ™ r 

4s = ,E d d *MM*)\ 3 + \vM*)\ a ) 

a=l 



+^S il ,..., ii ^(x)...^(x)} 

n . 

+ £ d d xd d yg(\x-y\)fc(x)% ) (y).(16) 

Here, the coupling of the replicas is given by the cor- 
relation function g(r) of Eq. ©, Greek indices denote 
replicas and the replica limit n — > is implied. 

For small k, the Fourier-transform g(k) of g(r) © 
reads: 

g(k) ~v + w \k\ a - d . (17) 

Thus, rewriting Eq. (|16fl in momentum space and tak- 
ing Eq. (|17|) into account, one obtains an effective La- 
grangean with three bare couplings uq, vq, wq. For a > d, 
the uio-vertex does not introduce additional divergences 
at k = and is irrelevant in the renormalization group 
sense [5(| H|j- The polymer limit m = leads to 
further simplifications. As pointed out in |4|, once the 
limit m, n — > has been taken, the uq and vq terms are 
of the same symmetry, and an effective Lagrangean with 
one coupling (uq + vq) of 0(mn = 0) symmetry <|13|) re- 
sults. This leads to the conclusion that weak quenched 
uncorrelated disorder i.e. the case a = d is irrelevant for 
polymers, and consequently also for star polymers. For 
a < d, the momentum-dependent coupling wok a ~ d has to 
be taken into account. Note that g(k) must be positively 
definite being the Fourier image of the correlation func- 
tion. This implies wq > for small k. Also, we assume 
the coupling uo to be positive, otherwise the pure system 
would undergo a 1st order transition. 

The resulting Lagrangean in momentum space then 
reads: 
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C ^ = 2 E E^o + * 2 )<£(*) + ^ E E 5 ( fcl + ■ • • + k *)$M • $M$ a (k 3 ) ■ $ a (k 4 ) 

a=l k ' a=l {k} 

n 

+ EE^ fcl + • • • +k4)\ki+k 2 \ a - d Mh) ■ Mk 2 )Mk 3 ) ■ $ {k 4 ). 

a/3 {k} 

Here, we have redefined uq + vq —> uq and denoted the scalar product by <fi ■ <fr. 
The replicated composite operator (|14f) reads in momentum space: 

n 

[<t>]f(ki, ...,kf) = S(k 1+ . . . +k f ) J2 Th,-,if€ ( fc i) • ■ ■ ^( fc /)- (18) 

a=l 



III. RENORMALIZATION GROUP APPROACH 

In order to extract the scaling behavior of the model 
(|18fl . and of the composite operator (|14|) we apply 
the field-theoretical renormalization group (RG) method 
[EE ■ We choose the massive field theory scheme 



with renormalization of the one-particle irreducible ver- 
tex functions T^' N '{ki, .., k^\pi, -,Pn\ MailAo}) a t non- 
zero mass and zero external momenta |76| . The one- 
particle irreducible (1PI) vertex function can be defined 
as: 



0(A) . . . ^R N ))^d d Ri . . . d d R N d d n . . . d d r L . 



(19) 



Here, {Ao} stands for the set of bare couplings uo,wq 
of the effective Lagrangean, {k}, {p} are the sets of ex- 
ternal momenta, Ao is the cutoff, and the averaging is 
performed with the corresponding effective Lagrangean, 
£ c g. To extract the anomalous dimensions of the com- 
posite operators l|14l) we define the additional /-point 
vertex function Ti , with a single [</>]{ insertion. Up to 
2nd loop order the graphs for r* can be derived from 
the usual graphs for I^ 0,4 ) by replacing in turn each four 
-point vertex by [</>]{ (see Fig. 3). 

The renormalized vertex functions T^' N ^ and are 
expressed in terms of the bare vertex functions as follows: 



r^ Ar) (W;M;A 2 ;{A}) = 
^< /2 r L ' w) ({fc};M;Ao;{A }), 
ri / i? ) ({fc};A 2 ;{A}) = 

Z„ / ri£ ) ({fc};/i§;{Ao}) ) 



(20) 



where Z^, Z^2, Z^f are the renormalizing factors, /t, {A} 
are the renormalized mass and couplings. 

The change of couplings uq^wq under renormalization 
defines a flow in parametric space, governed by corre- 



oo 



(/) 



FIG. 3: The graphs contributing to the vertex function F, 
up to 2-loop order, (a) represents the /-point vertex [<j>]l, (b): 
one-loop contribution, (c)-(f): two-loop contributions. 



sponding (3- functions: 
du 



dint 



, (3 w (u,w) 



dw 



dim 



(21) 



stands for evalu- 



where I is the rescaling factor, and 

ation at fixed bare parameters. The fixed points (FPs) 
u* , w* of the RG transformation are given by the solution 
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of the system of equations: 

f3 u (u*,w*)=0, 



p w (u*,w*)=0. 



(22) 



The stable FP, corresponding to the critical point of the 
system, is defined as the fixed point where the stability 
matrix 



B, 



(23) 



possesses eigenvalues {wi} with positive real parts. The 
flow of the renormalizing factors Z$, Z$a, Z*f in turn 
defines the corresponding RG functions: 



7<f>(u,w) 
7 2 (u,w) 

T)f(u,w) 



dlnZs 



d In Zm 



d\sxl 
d In Z^j 



dln£ 



(24) 
(25) 
(26) 



The critical exponents are the values of these functions 
Ig^l-© at tne stable accessible FP of Eq. (22J: 

v- X = 2-^(u*,w*)~%2(u*,w*), (27) 

r] = ^(u*,w*), (28) 

rj f = rif (u*,w*), (29) 

7/ = 1 + vr,j{vr,w*) + (^(2 -V)~ !)/• (30) 

Here, rjf is the anomalous dimension of the composite 
operator [</>]{. The expressions for the exponents v, r\ of 
a single polymer chain in long-range-correlated disorder 
we derived in Ref. 0,0, El- Only the RG functions rjf 
that correspond to the anomalous dimensions of the com- 
posite operator [(f>]l in the presence of correlated disorder 
remain to be calculated in order to extract the spectrum 
of star polymer exponents 7/ (given by Eq. I|3U|) ). 



IV. THE RESULTS 

The perturbative expansions for the functions l)2ip. 
(|24|l - (I26[) may be analyzed by two complementary ap- 
proaches: either by exploitin g a double expansion in 
e = 4-d,<5 = 4- a[ll HBl HI or by evaluating 
the theory for fixed values of the parameters d and a 
EE El Hiil. In the following we make use of both 
ways of analysis. 



A. One-loop approximation: e,8 - expansion 

For the qualitative analysis of the first order results, 

we apply a double expansion in e = 4 — d and 5 = 4 — a. 

First, we need to calculate the /-point vertex function 
(f) 

r* with a single insertion of the composite operator 



In the one-loop approximation we get: 



Tl J '(u ,w , {k} = 0,/io) = 1 - u 



dq 



(<7 2 + A§) 



2\2 



W ■ 



G 



fif - 1) 
6 

dqq a - d 



(31) 



We define renormalized mass fi 2 and couplings u, v by 
the renormalization conditions: 

A 2 = r^(fc,/i 2 ,M,u;)| fc= o, 

u = r i?,l(W)A 2 »'")^)l{fe}=o, 
w = r K,l({ fc }'A 2 ,«,^)l{fc}=o- 

The renormalization condition for the vertex function 
with [(f)] {-insertion is given by 



Z* f (u,w) = [rS^ («,«;)] -1 



= 1 + u 



m^) h _ w fif^i I: 



G 



G 



(32) 



Here, I\ and I2 are loop integrals given in the appendix. 
The expressions for the RG (3- and 7-functions (|21ll . 124|l . 
(|25|l within the same approximation read [T(ij |: 



(3 U = 





4 2r" 






£ 


U U Jl 


— S2uw 
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+ (25-e)-w 2 I 3 , 



W = -5 



w + — w I2 + w Di 



(33) 



+ e-[wuh], (34) 



1i IV IV 

14,2 =£ 3 Jl ~ S 3" /2 ' ^^ S ^ Dl - ( 35 ) 

Again, the loop integrals 1\ — /3,-Di are given in the 
Appendix. Note that contrary to the usual 4 theory the 
70 function in Eq. (|35() is nonzero already in the one-loop 
order. This is due to the fc-dependence of the integral D\. 
Combining Eqs. I|2G|) and l|21|) one defines r\j via familiar 
expressions l|3^|) and lp£f)l : 

??/ = Pu{u, w) d]l \ Z * f + P w {u,w) d ^ Z * f . (36) 



du 



dw 



To proceed with the analysis, we insert the expansions 
of the one-loop integrals: 



1 / £ 

h = - 1 1 - - 

e V 2 

1 



/ 2 



1 



1 



25 -£ 
1 /<5 
5 



28 — e 



(37) 
(38) 

(39) 

(40) 



Substituting (|37|l - l|4U|l into the expressions for (3- 
functions l|33|l , l|34|) and solving the FP equation (1221) , one 
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finds three fixed points: the Gaussian (u* = 0, w* = 0), 
the pure (u* — jS,w* = 0), and the non-trivial, long- 
range-correlated, LR fixed point: (u* = " "'" 



4 e-S ' 



W 



3 S(e-28)-. 



2 g_ e ). The analysis of the conditions of their sta- 
bility and accessibility we performed in Ref. jlflj . The 
results are displayed schematically in Fig. at e, <5 > 0, 
the crossover from the pure FP to the LR takes place at 
8 = e/2, i.e. a = 2 + d/2. Note, however, that the LR FP 
is stable in the region a > d, where the influence of the 
disorder is expected to be irrelevant, see the explanation 
after Eq. (|17|l . These first order results give a qualitative 
description of the crossover to the new universality class 
in the presence of long-range-correlated disorder. 

The expression for the critical exponent v reads |10| : 



i/(Q) = l/2 + e/16, 
z/W = 1/2 + 5/8, 



8 < e/2, 
e/2 < 8 < e. 



From (|36f) we find: 




|e/(/-l), <5<e/2, 
iSf(f-l), e/2<S<e, 



(41) 



(42) 



the star exponent 7/ at fixed / > 3, when the correla- 
tion of the disorder becomes stronger (i.e. parameter a 
decreases). However, the behavior for chain polymers i.e. 
for / = 1,2 differs: in this case the exponents 71 = 72 
increase for decreasing a. 




f 



FIG. 5: Comparison of the exponents 7/ for star polymers 
in a pure solution (continuous line) and in a solution inside 
a correlated porous medium (broken line), following Eq. l43H 
for e = 1, 8 = 0.9. 



W 



U 



(a) 



w 



u 



(b) 



FIG. 4: Fixed point picture for d < 4, a > d. Stable physically 
accessible fixed points are shown by squares the unstable ones 
by discs, (a): 8 < e/2, the pure fixed point (u* / 0,w* — 0) 
is stable. At S = e/2, it interchanges its stability with the 
LR fixed point {u* / 0,w* 0). (b): for e/2 < 8 < e the 
LR fixed point becomes physically accessible and stable. The 
Gaussian fixed point (u* — 0, w* — 0) is stable for d > 4, a > 
4. 

Substituting (IT) into pD)l. finally we get: 



7/ 




1 w - 



-3) 
3) 



6 < e/2, 
e/2 < 8 < e, 



(43) 



The first line in l|43() recovers the exponent for the /-arm 
polymer star in the pure solution |3j| , whereas the second 
line brings about a new scaling law. 

To obtain a naive estimate of the numerical values of 
these exponents, one can directly substitute into Ij43(l the 
value e = 1 (corresponding to d — 3) and different fixed 
values for correlation parameter a We note a decrease of 



This crossover is also clearly seen in Fig. [S] where we 
compare the behavior of 7/ for the case with and without 
correlated disorder. As this figure shows, the correlation 
of the disorder effectively enhances the deviation from the 
mean field value 7^ = 1 which is positive for / = 1 , 2 
and negative for / > 3 in this approximation. 



B. Two-loop approximation: fixed d, a approach 

To obtain a quantitative description of the scaling be- 
havior of star polymers in long-range-correlated disorder, 
we proceed to higher order approximations. We make 
use of the fixed d — 3 RG approach jT^ , considering the 
massive RG functions at fixed space dimension d. Also 
the additional parameter a in the expansions for the RG 
functions in renormalized couplings u, w (|21|l . I|24|l - H26() 
is fixed in this approach and we work hereafter with these 
expansions. As is well known [liJ 03, j such expan- 
sions are in general characterized by a factorial growth 
of the coefficients which implies a zero radius of conver- 
gence [73 ■ No reliable data can be extracted from a naive 
analysis. For the present model, this particular feature 
shows up already in the first order of perturbation theory 
in u and w. Indeed, for the plain one- loop f3- functions 
(|34|l a non-trivial FP LR does not appear if one solves 
the non-linear fixed point equation 1)22(1 directly at d = 3 
and 2 < a < 3. To take into account higher order contri- 
butions, the standard tools of asymptotic series resum- 
mation have to be applied fn\ . 

The two-variable Pade-Borel resummation technique 
|78| | that we use consists of several steps. Consider the 
two-variable series for a RG function h(u,w). First, we 
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construct the Borel image of the initial function: 



h(u, w) = ^2 at 



E 



r'(i + j + i) ' 



where T(i + j + 1) is Euler's gamma function. Then, the 
Borel image is extrapolated by a rational Pade approxi- 
mant [K/ L](u,w). This ratio of two polynomials of order 
K and L is constructed as to match its truncated Tay- 
lor expansion to that of the Borel image of the function 
h{u, w). The resummed function is then recovered by an 
inverse Borel transform of this approximant: 



h TCS (u,w)= / dtcxp(-t)[K/L](ut,wt). 



(44) 



In our previous work [l2L Il3| we have analyzed the 
resummed expressions for the two-loop RG functions of 
the model of a single polymer chain in the long-range- 
correlated disorder in three dimensions, and found that 
a fixed point LR appears and is stable at 2.2 < a < 3. 
This FP disappears at a < 2.2 and the pure SAW FP 
remains unstable. This behavior may be interpreted to 
indicate, that the presence of stronger correlated disorder 
(at a < 2.2) might lead to a collapse of the polymer chain. 
To obtain a quantitative picture of the scaling behavior 
of star polymers, we only need to extend these results 
by a calculation of the renormalization factor Z*f 
Taking into account the two-loop contributions shown in 
Fig. |2| we get: 



Z 



*f 



l + u h - w 1 2 

6 6 

W /(/ - 1) \(-^(f - 2)(/ - 3) - 1) + \ ) /f 



/ " 2 + 1 I / 



-uwf(f-l) 



~(/-2)(/- 3) + !/(/-!) h 



Di 



+™ 2 f(f - 1) 



(/-2)(/-3) f(f-l) 



72 



/-2 l\. f-2 1 1 2 

— + 6) h ~ + 9 hh 18 h 3/(/-l) 



(45) 



r 



The expressions for the loop integrals I\ , . . . , iio , D\ and 
their numerical values at d = 3 and different a are pre- 
sented in the Appendix. 

In this way, the function rjf can be found, using l|3()|) 



and familiar expressions for the two-loop /3-functions as 
given in 03, |l3( . The resulting two- loop expansion for 



J 



/(/-l) (4-a)/(/-l) 7 
-u w h/h 



16 



32 32 



(46) 



I I -^(l + (4-a)) + ^(l-(4- fl )) + (4- a )^/ 1 / 2 
+ (( 1 + ( 4 - a ^S-¥ + ^)) (/7 + /9 ) 



16 16 Jf(f-l)' 

, (4-o)/(/-l) 



16 



-Is 



9 



, (4-g)(/-l)(/ + l) T , ( 4 -g)/(/-l)(/-2) 
+ 16 h + 16 



Inserting the series for v and r\ for the polymer chain 
in the long-range-correlated disorder from Refs. [l2l 
together with r/f of Eq. I)46|l into (|3UII . we finally obtain 
the corresponding scries for 7/. Substituting the numer- 
ical values of the LR correlated FP, found for different a 
from Refs. 0, ^| and applying a Pade-Borel resumma- 
tion as explained above we get the numerical values for 
exponents 7/ in three dimensions for different values of 



the correlation parameter a and number of arms /. Our 
final estimates that result from this procedure are pre- 
sented in Table For / — 3, the first order contribution 
to 7/ is zero, whereas it is non-zero for 773. Therefore, 
to obtain a resummed value for 73 we have resummed 
the series for 773 using the values for r\ and v of chain 
polymers in long-range-correlated disorder. 



a\f 


i;2 


3 


4 


5 


6 


7 


8 


9 


3, [42, 43] 


1.18 


1.06 


0.86 


0.61 


0.32 


-0.02 


-0.4 


-0.8 


3, [36] 


1.1573(2) 1.0426(7) 0.8355(10) 0.5440(12) 0.1801(20) 


-0.2520(25) 


-0.748(3) 


-1.306(5) 


3 


1.17 


0.99 


0.83 


0.57 


0.26 


-0.08 


-0.56 


-0.87 


2.9 


1.25 


0.87 


0.78 


0.46 


0.09 


-0.32 


-0.76 


-1.23 


2.8 


1.26 


0.81 


0.76 


0.43 


0.06 


-0.36 


-0.80 


-1.26 


2.7 


1.28 


0.74 


0.72 


0.40 


0.01 


-0.40 


-0.85 


-1.31 


2.6 


1.30 


0.73 


0.70 


0.37 


-0.03 


-0.46 


-0.91 


-1.37 


2.5 


1.34 


0.71 


0.70 


0.35 


-0.10 


-0.51 


-1.00 


-1.44 


2.4 


1.35 


0.70 


0.70 


0.31 


-0.10 


-0.55 


-1.02 


-1.50 


2.3 


1.38 


0.70 


0.69 


0.29 


-0.13 


-0.59 


-1.06 


-1.55 



TABLE I: Critical exponents 7/ for the /-armed star in three dimensions at different values of the correlation parameter a. 
The first and the second rows (a = d = 3) present results for a polymer star in a good solvent without porous medium obtained 
within the field-theoretical RG in three- loop approximation, Ref. |42l I4H . and by the Monte Carlo simulations, Ref. |36jl . 
correspondingly. 



Let us recall, that for a = d = 3 the problem is equiva- 
lent to the situation without structural disorder. There- 
fore, in the first two rows of Table[I]we give RG estimates 
for the exponents 7 f obtained in a three-loop approxima- 
tion in Ref. |42l Eflj as well as recent data of MC sim- 
ulations [3^ | . Comparing these data with our two- loop 
results (the third row of the Table) allows to estimate 
the consistency of the calculational scheme that we ap- 
ply. The good mutual agreement found for the low values 
of / supports our approach. The fact that the discrep- 
ancy increases with / is expected, taking into account the 
strong combinatorial /-dependence of the coefficients of 
expansions 145(1 . (|46(1 . This growth is difficult to control 
in a consistent way during the resummation. 

As we noted above, the choice / = 1, / = 2 recovers the 
case of a single polymer chain. Therefore, the first and 
the second columns of Table 1 give an estimate for the 
dependence of the exponent 7, Eq. (1): = 7^ = 



72 . The remarkable feature of the estimates for 7^ 
listed in the Table 1 is that they predict a qualitatively 
different behavior of jy' for / = 1, / = 2 and / > 3. 
Indeed, as one sees from Table 1, a decrease in a leads 
to an increase of J^,^^ while jf for / > 3 decrease 
in this case. This tendency is also found for the one-loop 
e, (^-expansion (Eq. JSJl). 

Recall, that the scaling exponent of a star polymer in 
a pure solvent is given by 7| a_3 ^ and let us return back to 
Eqs. J7J) and JSJ for the free energy of a star in the pure 
solvent and in a porous medium. Then our results indi- 
cate two different regimes of the entropy-induced change 
of the polymer concentration for a solvent in a porous 
medium with respect to the pure one. Namely, the free 
energy of the chain polymers (/ = 1,/ — 2) is reduced 
by the presence of correlation in a porous medium. On 
the other hand, the free energy of a star polymer (/ > 3) 



is increased by correlations of the environment. 
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2.9 
2.8 
2.7 
2.6 
2.5 
2.4 
2.3 



641 643 044 04 



e 



40 



e 4 



e 4 



e 4 



1.306 1.627 1.865 2.042 2.220 2.348 2.402 2.531 

1.286 1.565 1.777 1.932 2.071 2.163 2.246 2.345 

1.262 1.502 1.691 1.817 1.941 2.024 2.082 2.166 

1.239 1.459 1.608 1.739 1.843 1.910 1.987 2.040 

1.229 1.410 1.554 1.668 1.762 1.834 1.876 1.929 

1.217 1.392 1.521 1.617 1.705 1.758 1.799 1.883 

1.193 1.360 1.474 1.574 1.651 1.707 1.725 1.772 



TABLE II: Contact exponents 0//, governing the scaling be- 
havior of interaction force between two /- and /'-armed poly- 
mer star in three dimensions at different values of the corre- 
lation parameter a. 




FIG. 7: The contact exponent 0//' as function of / and /' 
for fixed correlation parameter a = 2.7atd = 3. 



To investigate the influence of a porous medium on 
the effective interactions between star polymers we cal- 
culate the contact exponents 9//' (Eq. UJ. Our results, 
obtained by a Pade-Borel resummation of the series de- 
rived from Eq. 10. are presented in Figs. El Hand for a 
selected set of exponents also in Table [H] In Fig. 5 we 
show the contact exponent Off for two stars of the same 
number of arms / as a function of a and /. The exponent 
increases with increasing of / and a for / > 3. Fig. 6 
presents O//' for a fixed a (we have chosen a = 2.7 for 
an illustration). For fixed /, /', this exponent decreases 
with the decrease of the correlation parameter a. Thus, 
we can conclude, that polymer stars interact more weakly 
in media with strong correlated disorder. 




FIG. 6: The contact exponent Qff as function of / and corre- 
lation parameter a at d — 3. Each line shows the dependence 
of 0// on a at fixed /. 



V. CONCLUSIONS 

The present study provides numerical estimates for the 
spectrum of critical exponents that govern the scaling be- 
havior of the /-arm star polymers in a good solvent in 
the presence of a correlated disordered medium, charac- 
terized by a correlation function g(r) ~ r~ a at large dis- 
tances r. This extends previous results GHEUTa that 
have shown that the scaling behavior of polymer chains in 
this type of disorder belongs to a new universality class. 

Working within the field-theoretical RG approach, we 
applied both a double expansion in e = 4— d and 5 = 4 — a 
as well as a technique that evaluates the perturbation se- 
ries for fixed d, a. The first one-loop analysis allowed us 
to identify a quantitatively new behavior in comparison 
with the pure case. The second approach, refined by a 
resummation of the resulting divergent series, resulted 
in numerical quantitative estimates for the scaling expo- 
nents. We found the numerical values of the exponents 7/ 
in three-dimensional case for different fixed values of the 
correlation parameter 2.3 < a < 2.9, and for fixed num- 
bers of arms / = 1, . . . , 9. Depending on the value of 
/, we find two different regimes of the entropy-induced 
effects on the polymer in a correlated porous medium. 
While an increase of the correlation of the disorder causes 
the free energy of chain polymers (/ = 1,/ = 2) to de- 
crease, the same change in correlation rather leads to an 
increase in the free energy for star polymers (/ > 3). 
Therefore, for a mixture of chain and star polymers of 
equal molecular mass (same total number of effective 
monomers) in a solution for which a part of the volume 
is influenced by a porous medium the disorder-influenced 
part of the solvent is predicted to be enriched by chain 
polymers. Correspondingly, the relative concentration of 
star polymers to and chain polymers will be lower in the 
porous medium. 
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From our numerical estimates for contact exponents 
Off, we deduce the influence of the correlated disorder 
for the effective interaction between star polymers. Again 
we find different behavior for chain and star polymers. 
While for chain polymers the effective contact interac- 
tion increases for decreasing a, i.e. for enhanced corre- 
lation, the mutual interaction between star polymers is 
weakened in correlated media. 
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VI. APPENDIX 

Here, we present the expressions for the loop inte- 
grals, as they appear in the RG functions. We make 
the couplings dimensionless by redefining u — ujj, d ~ 4 and 
w = u>/i a ~ 4 . Therefore, the loop integrals do not explic- 
itly contain the mass. 



{q 2 + l) 2 



h = 



h = 



h = 



h = 



'10 



dqq a ~ d 
(g 2 + l) 2; 

d[j - q 2(a-d) 
(q 2 + l) 2 1 



dqidq 2 q[ 



(a-d) 



(d + l) a ((9i-©) 2 + l) s 
dq 1 dq 2 q[ a -Vq a 2 - d 



(<zf- 


hl)2((g 1 - ffiJ )2 + 


l) 2 ' 






dq\dq 2 






(ql- 


-l)( 92 2 + l) 2 (( gi - 


- q2? - 






dqidq 2 q?- d 






{ql- 


h l)( g | + l)2(( gi . 


- q2? - 


HI)' 




j — ♦ j — ♦ a—d (a 

dqidq 2 q x q\ 


-d) 




(qf- 


h l)( 5 | + l)2(( gi . 


- q2? - 


hi)' 




dq x dq 2 q\~ d 






(9?- 


~l) 2 (q 2 + l)(( qi - 


- q2? - 


-1)' 




dqidq 2 q{- a ~ 


d) 





-i) 2 fe 2 + i)((gi-g 2 ) 2 + i)' 

dq q a 



2 

-d 



[g + fc]2 + l)j fc2=0 - (4?) 
The correspondence of the integrals to the diagrams in 
Fig. El is: (b): integrals h, h, h, (c): h,h; (d) : hh, 
7|; e: I 6 , J 7 , I 8 , h, ho; (/): h,h,h,h. In our calcu- 
lations, we use the following formulas for folding many 
denominators into one (see e.g. [58]): 



£(ai+ . . .+«„) 
r(ai) . . .T(a n ) 



dx\. . . I dx n 



-xi- 



-X n -l) 



{x\a\+ . . . +a;„_ia n _i+(l— x\- 



-x n -i)a n ) ai 



«l+..-+o„ 



(48) 



To compute the d-dimensional integrals we apply 

oo 

dqq*- 1 



2\ot 



(q 2 + 2kq + m 2 ) 

lmma-d/2) /2 _ a 

2 r(a) 

As an example we present the calculation of the inte- 
gral I 7 . First, we make use of formula 148(1 to rewrite: 



1 



(<? 2 + !)((<?! -<Z 2 ) 2 + 1) 



da; 



r(2) 

r(l)r(l) J (q 2 + 2xqiq 2 + xq 2 + l) 2 ' 



(50) 



Now one can perform integration over q\ , passing to the 
d-dimensional polar coordinates and making use of the 
formula (|4*5|) : 



dqi ft 



a—d 



(qf + 2xq[q 2 



xq\ 



C 



dqi qf 



(q 2 + 2xq x q 2 + xq 2 + l) 2 
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a 

2.9 
2.8 
2.7 
2.6 
2.5 
2.4 
2.3 



h 
0.7855 
0.7855 
0.7855 
0.7855 
0.7855 
0.7855 
0.7855 



h 
0.7155 
0.6605 
0.6170 
0.5825 
0.5550 
0.5345 
0.5185 



h 
0.6605 
0.5825 
0.5345 
0.5085 
0.5000 
0.5085 
0.5345 



h 
0.5621 
0.5187 
0.4848 
0.4575 
0.4361 
0.4198 
0.4075 



h 
0.5119 
0.4363 
0.3807 
0.3393 
0.3080 
0.2857 
0.2688 



h 
0.4114 
0.4114 
0.4114 
0.4114 
0.4114 
0.4114 
0.4114 



h 
0.3643 
0.3274 
0.2981 
0.2746 
0.2555 
0.2406 
0.2283 



h 
0.3477 
0.3016 
0.2677 
0.2425 
0.2238 
0.2106 
0.2014 



h 
0.3916 
0.3756 
0.3626 
0.3525 
0.3448 
0.3395 
0.3365 



ho 
0.3756 
0.3525 
0.3395 
0.3357 
0.3408 
0.3557 
0.3823 



Pi 
0.0052 
0.0012 
0.0015 
0.0021 
0.0027 
0.0034 
0.0041 



TABLE III: The numerical values for loop integrals Ii for d = 3 and different values of the correlation parameter a. 



j gyS (| + ^., )p M, m 

where the constant C = (27r) d / 2 /r(d/2) results fr om in- 
tegration over the angular variables. It does not appear 
explicitly in the following expressions. Finally, we are 
left with: 

h = ir(a/2)r(2-o/2) 



r i mi + qix{i " x))a/2 " m 



this integral was calculated numerically, fixing the values 
of the parameters d, a using the MAPLE package. Note 
that some of the integrals can also be evaluated analyti- 
cally. 



Below, we list the results for all the integrals. 



h = -r(d/2)r(2 - d/2), 

h = ir(a/2)r(2-a/2), 

h = ir((2a-d)/2)r(2-(2a-d)/2), 

h = ^r(a/2)r(d/2)r(2-a/2)r(2-d/2), 

7 5 = -r(a/2)r(o/2)r(2-o/2)r(2-a/2), 

h = ir(d/2)r(2 - d/2) x 

oo j d—1 f 



7 8 - -r(a/2)r(2 - a/2) x 



o + 



^ y dx(l + g 2 2 x(l - i)) 



a/2-2 



7 9 = -r(d/2)r(2 - d/2) x 



o (9 2 2 + l) 2 



&x{l + q 2 2 x{l - x)) d ' 2 - 2 , 



1, 



Jio = -r(d/2)r(2 - d/2) 

^ 1 



oo i _2a— d— 1 



o 

(a-2)(a-3)(a-4) 
192sin(7r(l/2a - 1)) ' 



Note, that the analytical value for D\ is taken from [67 
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